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ABSTRACT 

We report the discovery of two low-luminosity quasars at z ~ 4, both of which 
show prominent N IV] A1486A emission. This line is extremely rare in quasar 
spectra at any redshift; detecting it in two out of a sample of 23 objects (i.e., 
~ 9% of the sample) is intriguing and is likely due to the low- luminosity, high- 
redshift quasar sample we are studying. This is still a poorly explored regime, 
where contributions from associated, early starbursts may be significant. One 
interpretation of this line posits photoionization by very massive young stars. 
Seeing N IV] A1486A emission in a high-redshift quasar may thus be understood 
in the context of co-formation and early co-evolution of galaxies and their su- 
permassive black holes. Alternatively, we may be seeing a phenomenon related 
to the early evolution of quasar broad emission line regions. The non-detection 
(and possibly even broad absorption) of N V A1240A line in the spectrum of 
one of these quasars may support that interpretation. These two objects may 
signal a new faint quasar population or an early AGN evolutionary stage at high 
redshifts. 

Subject headings: quasars: emission lines — quasars: general — galaxies: evo- 
lution — quasars: individual (DLS1053— 0528) — quasars: individual (ND- 
WFS1433+3408) 
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Introduction 



A growing body of theoretical models backed by emerging observational evidence paints 
a picture of joint formatio n a n d evo l ution of gal a xies a nd quasars; for reviews see, e.g., 
the proceedings edited by iHol (120041 ). iDjorgovskil (120051 ). and references therein. In this 
picture, merger-driven build-up of massive ha los trigger both vigorous star formation and 



fuel the central black hole igniting a quasar (IKauffmann fc Haehneltl l2000l : iHopkins et al. 



20061 ). AGN feedback is now believed to be an essential factor in the formation and evolution 
of galaxies. Understanding the AGN-starburst connection is especially interesting at high 
redshifts, as we probe the epoch of the initial assembly of massive galaxies. Observationally, 
this is challenging in the case of host galaxies of luminous Type I quasars, as the AGN 
activity far outshines any star formation related processes. However, by observing lower 
luminosity Type I quasars or Type II (e.g., obscured) quasars, the blinding intensity of the 
central engine is reduced, allowing properties of the host galaxy to be studied. 

In this Letter, we report the discovery of two low-luminosity Type I quasars where 
we may be seeing evidence for a simultaneous starburst with a top-heavy IMF, and AGN 
activity. Both quasars show a prominent, moderately-broad N IV] A1486A emission line, 
which is rarely seen in quasar spectra at any redshift. 



2. Observations 



The two objects are low-luminosity quasars at 2; ~ 4, found in a sample of 23 objects 
used to measure the low-l uminosity end of the quasar luminosity function (QLF) at 2; ~ 4 
( iBogosavljevic et al.l 120071 ). The parent sample ranges in redshift from z = 3.7 to z = 5.1. 
The candidates for this sample were selected from the NO AO Deep Wide Fie ld Survey (^ND 



WFS) Bootes field (jjannuzi fc Deylll999l ) and the Deep Lens Survey (DLS: IWittman et al. 
2OO2I ). Quasar candidates were selected based on the colors of simulated quasars in the 
3.5 < z < 5.2 redshift range, incorporating a range of spectral slopes, Lya equivalent 
widths, and intervening neutral hydrogen absorbers, down to the limiting magnitude of 
R = 24. Finding charts for the two N IV] quasars are presented in Figure [H Details of the 
survey and candidate selection will be presented by Bogosavljevic et al. (2007). 

We obtained spectroscopic followup for our candidates on UT 2006 M ay 20 through 
22 with the Low- Resolution Imaging Spectrometer (LRIS; lOke et al.l Il995l ) on the Keck 
I telescope. Only the red camera was used, with the 400 lines mm~^ grating blazed at 
8500A. The spectra were reduced using BOGUS, an IRAF package developed by Stern, 
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Bunker, & StanforcO for reducing slitmask data, modified slightly for our single-slit data. 
DLS1053— 0528 was discovered in a 1800 second exposure of a candidate in the DLS F4 
field and NDWFS 1433+3408 was identified from a 900 second exposure of a candidate in 
the Bootes field. The final spectra are presented in Figure [2] and the source parameters are 
listed in Table [H 

Several things are worth noting. First is the detection of the N IV] A1486A emission 
lines, which is particularly strong in DLS1053— 0528 (e.g., in comparison with the C IV 
A1549A line). At the same time, the commonly observed N V A1240A emission line is 
entirely absent, and may even be seen in absorption in DLS1053— 0528 (see Figure [3]). This 
latter observation is somewhat unusual since the C IV A1549A line in this source does not 
show broad absorption. While broad absorption lines are seen in approximately 10% of 
quasars, such quasars will show broad absorption in all permitted species. 

The permitted emission lines for the two N IV] quasars are narrow compared with 
typical quasar line-widths. A single component Voigt fit to Lya measures a full-width at 
half-maximum (FWHM) of ~ 500 km s~^. Since the blue side of the Lya profile is absorbed, 
we forced symmetry in the line by mirroring the red side of the line profile over the peak 
wavelength. To search for a broad-line component we fit a two-component Gaussian to 
the Lya profile. For DLS1053— 0528, we measure a narrow-line FWHM of 434 km s~^ 
and a broad-line FWHM of 1727 km s~^. In NDWFS 1433-1-3408 we measure a narrow- 
line FWHM of 713 km s~^ and a broad-line FWHM of 3015 km s~^. Therefore, while 
these objects have significant contribution from their narrow-line components (50% of the 
line flux in DLS 1053-0528 and 24% in NDWFS 1433+3408), their broad-line velocities of 
> 1000 km s~^ place them in the quasar regime. Using the velocity widths from the broad- 
line component of Lyg a s a pro xy fo r C IV or H/j, and the black hole mass estimators 



from iDietrich fc HamannI (120041 ) and IVestergaard fc Peterson! (120061 ). we obtain Mbh 



(5-15) X 10*5 Mq for DLS1053-0528 and Mbh ~ (32-77) x 10^ Mq for NDWFS1433+3408. 
Consistent with the somewhat narrow line widths (for a broad-lined quasar) and the low- 
luminosity, these inferred black hole masses are much lower than the masses inferred for 
luminous quasars. 

On UT 2007 January 23 we imaged DLS1053— 0528 with the Slit-viewing Camera 
(SCAM) of the NIRSPEC instrument on Keck-II. We used the J-band filter and imaged 
the source with a 9-point dither pattern of 120 second exposures, for a total of 1080 sec- 
onds at the position of the quasar. The seeing was 0'.'7. The data were reduced using the 
XDIMSUM package in IRAF. We calibrate our image to the 2MASS J-band using a 15.26 



^ https : //zwolf kinder . jpl . nasa . gov/$\sim$stern/h.omepage/bogus . html 



-4- 



magnitude point source that is detected in the field. The quasar is not detected, with a 3a 
magnitude threshold of 24.9 (Vega) magnitudes. 

NDWFS1433+3408 was obs erved as part of F LAMEX, a near-infrared survey overlap- 
ping 4.1 deg^ of the Bootes area (lElston et al.ll2006l ). The quasar is detected at ~ 5a in both 
the J and Kg images. In addition, this quasar was ob served in the IRAC Sh allow Survey, 
an 8.5 deg^ Spitzer imaging survey of the Bootes field (lEisenhardt et al.ll2004l ). The source 
is faintly (< 5a) detected at 3.6/im and 4.5/im at 19.0 and 19.1 (Vega) magnitudes, respec- 
tively. We list the 4" aperture magnitude from the FLAMEX catalog as well as the IRAC 
detections in Tabled! The quasar is not detected in the 5.8/xm or 8.0/im images, whose 5a 
detection thresholds are 15.9 and 15.2 (Vega; 3" aperture) magnitudes, respectively. Note 
that the [3.6] — [4.5] colors for this ob ject are quite blue compared to typical mid-infrared 
colors of AGN (e.g.. lStern et al.ll2005l ). The result is not unexpected: at z ~ 4, Ha emission 
is shifted into the shortest wavelength IRAC band, causing AGN at this redshift to have 
blue [3.6] — [4.5] colors. 

Figure m plots the spectral energy distributions (SEDs) of both object s, using all avail- 
able p hotometry. We compare this SED with the SPSS composite spectrum (jVanden Berk et al 
200ll ) as well as the starburst template from lKinney et al.l (119961 ) with E{B — V) < 0.1, both 
shifted to z = 3.88. Significant deviations are seen, especially at long wavelengths where 
both objects appear much bluer than the average quasar and starbusrt spectra. The SED 
of NDWFS1433+3408 seems more consistent with the starburst spectrum, at least in the 
rest-frame ultraviolet, but DLS1053— 0528 deviates strongly from both. 



Discussion 



The most striking feature in the spectra of these two faint quasars is the presence of 
the extremely unusual N IV] A1486A emission line. The quasar population generally shows 
remarkable spectroscop ic similar ity out to the highest observed redshifts, with no obvious 
signs of evolution (e.g.. lFanll2006l ). None of the published average quasar spectral templates 



show any trace of the N IV] A1486A line (Figure [S]). The top spectrum (solid line) is the 
SDSS quasar composite from lVanden Berk et al.l (120011 ). The quasars that contribute to this 
part of the spectrum are at redshifts comparable to the quasars in our sample, but they 
sample the bright end of the quasar lumin osity function . The bottom spectrum (dotted line) 
is the HST UV composite spectrum from iTelfer et al.l (|2002| ). The objects contributing to 
this part of the spectrum are low-redshift quasars {z < 1). 



The only strong detection of this line in a quasar that we are aware of is in the nitrogen 
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Tuasar Q0353— 383 (lOsmer fc SmithlllQSd ). This object has been analyzed by [Baldwin et al. 
[ I2OO3I ) who derive an overabundance of nitrogen by a factor as high as 5 — 15 times solar 



i n thi s object. A sample of apparently nitrogen-rich quasars was compiled by iBentz et al. 



( 120041 ) from the 6650 quasars in the SDSS DRl database with 1.6 < z < A.l (allowing for the 
detection of N IV] A1486A and N III] A1750A). Isentz et all kooi \ estimate that "nitrogen- 
enriched" quasars make up at most 0.2% — 0.7% of quasars. Their sample has luminosities 
in the range Mj = —28.07 to Mj = —24.63. We found two such objects in a sample of 23 
quasars (8.7%) that are concentrated at 2; ~ 4 and are 0.7 magnitudes deeper than the SDSS 
sample. They found 20 quasars with equivalent widths > 2.0A in both nitrogen lines and 
33 quasars with EWs > 2.0A in only one of the nitrogen lines. In these quasars the N IV] 
A1486A line is accompanied by similarly strong N V A1240A and N III] A1750A, the latter 
of which is outside our wavelengths range. We fit Voigt profiles to the nitrogen lines in our 
spectra (as well as to Lya and C IV A1549A) to determine line fluxes, equivalent widths 
and dispersion velociti es (Table [21 Figure [3]). The mean N IV] A1486A equivalent widths for 
the iBentz et al.l ( 20041 ) sample is 3.7A {aEw = 1.5A), while the equivalent widths of N IV] 
A1486A for our quasars are 28O.2A for the DLS source and 24. 5A for the Bootes source. 
This, combined with the non-detection (and possible absorption) of N V A1240A in the DLS 
source suggests that these quasars may be of a different ilk than the "nitrogen-enriched" 
population. 

In galaxies with no AGN signature, we are aware of only two other instances of this line 
being seen in the high-redshift Universe. One especially telling d etection is in the spe ctrum 
of the Lynx arc, a gravitationally lensed HII galaxy at z=3.357 (iFosbury et al.ll2003l ). The 
line intensities in this object's spectrum, strong N IV] AA1483, 1487A, O III] AA1661, 1666A, 
and C III] AA1907, 1909A, as well as the absence o f N V A1240, f avor a hot (80,000 K) 
blackbody over an AGN as the ionizing source. The iFosbury et al.l (120031 ) modeling of the 
spectrum suggests a top heavy IMF, which is especially interesting since such an IMF is now 
believed to be characteristic of ear ly, metal poor star formation, e.g. Population III stars. 
Alternatively, iBinette et al.l (120031 ) suggest an obscured AGN as the photoionizing source 
of the Lynx arc. Their model invokes dense absorbing ga s near a centra l AGN that filters 
the powerlaw and mimicks the hot blackbody suggested by IFosbury et al.l (120031 ) This model 
does require at least a weak N V A1240A line detection. 

The second known N IV] A1486A emi tter is a 2: = 5.55 galaxy in the GOODS survey 
( IVanzella et al.ll2006l : iFontanot et al.l 120071 ). Similar to the Lynx arc, this object shows an 
extremely strong Lya line as well as N IV] A1486A, while N V A1240A is absent. No detailed 
analysis has been published on this source as yet. 



An alternative possibility is that we are witnessing an evolutionary effect in the quasar 
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broad-emission-line region. This is suggested by the traditional broad-hne shape of N IV 
A1486A in DLS1053— 0528 as well the peculiar absence (or even absorption) of N V A1240A 
in its spectrum. Detailed modeling of such effects is beyond the scope of this paper. 

The detection of such a rare emission line in two out of 23 2; ~ 4 quasars in our sample 
(i.e., ~ 9% of the sample) suggests that it occurs more commonly in low luminosity quasars 
at high redshifts, a regime which we are now exploring systematically for the first time. At 
these luminosities we are probing deep into the QLF where the effects of a luminous star- 
bursts can be detected more easily; if there is an evolutionary trend, such phenomena may 
be present at high redshifts, and not in the better studied quasar samples at z ~ — 2. 
As the sample sizes of comparably faint AGN at these redshifts increase, we will be able to 
determine whether this is indeed a systematic evolutionary effect related to the early stages 
of co-formation and co-evolution of galaxies and AGN. 



We thank the referee for helpful feedback. We are grateful to the staff of W. M. Keck 
observatory for their assistance during our observing runs. This work was supported in part 
by the NSF grant AST-0407448, and by the Ajax foundation. The work of DS was carried 
out at Jet Propulsion Laboratory, California Institute of Technology, under a contract with 
NASA. 
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Table 1. N IV] Emitting Quasars 





R.A. 


Dec. 




R 


J 




[3.6] 


[4.5] 


Name 


(J2000) 


(J2000) 


z 


(mag)'^ 


(mag) 


(mag) 


(mag) 


(mag) 


DLS1053--0528 


10 53 46.1 


-05 28 59 


4.02 


23.83 


>26.4 








NDWFS1433+3408 


14 33 24.5 


+34 08 41 


3.88 


22.62 


21.49 


19.94 


19.0 


19.1 



"^AB magnitudes. All other magnitudes are Vega-based magnitudes. 
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Table 2. Quasar Line Properties 



Line Flux rest-frame EW FWHM 

Line (10~^^ erg cm~^ s~^) (A) (km s~^) 

DLS1053-0528: 

Lja 2.92 ±0.08 762.1 1726.8^ 
NV 

NIV] 1.048 ±0.008 280.2 807.5 

CIV 0.55 ±0.02 91.1 1101.5 

NDWFS1433+3408: 

Lya 3.8 ±0.3 244.6 3015.6^ 

NV 0.25 ±0.18 11.6 1086.9 

NIV] 0.21 ±0.07 24.5 524.9 

CIV 1.3 ±0.3 91.1 1334.0 



^These line widths are from the broad component of a 
two-component Gaussian fit to the Lya line. See §2 for 
details on the fitting procedure. 
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Fig. 1. — 3'x 3' finding charts from the i?-band images of the original survey data. 
DLS1053— 0528 is on the top (the gray stipe through the center of the image masks a bleed 
trail from a nearby saturated star) and NDWFS1433+3408 is displayed on the bottom. In 
these images, north is up and east is to the left. The white circles mark the locations of the 
quasars. 
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Ly« DLS1 053-0528 z = 4.02 




6000 6500 7000 7500 8000 

Wavelength (A) 

Fig. 2. — Keck LRIS spectra of the N IV] A1486A emitting quasars. Top: z = 4.02 quasar 
found in the DLS field. This spectrum lacks N V A1240A in emission but the hue appears 
present as a broad absorption hue. Bottom: z = 3.88 quasar found in the NDWFS Bootes 
field. This object shows a weaker N IV] A1486A line and a possible weak detection of N V 
A1240A. 
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NV NIV] 
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Fig. 3. — Voigt profile fits (dotted lines, with continuum model shown by the dot-dashed 
hne) to Nitrogen hues in our quasars. The hne fits to DLS1053— 0528 are shown in the top 
panels and the fits to NDWFS 1433+3408 are shown in the bottom panels. 
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Fig. 4. — Rest-frame spectral energy distribution (SED) for NDWFS 1433+3408 (squares) 
and DLS1053— 052 8 (triangle s) com pared with the SDSS quasar composite spectrum from 
Vanden Berk et ahl (so lid line: \200l\) and a n unobscured starburst template with E{B — V) < 
0.1 from lKinney et al. (dashed line; 1996). 
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Fig. 5. — Composite quasar spectra showing the a bsence of N IV] A 1486 A. Th e top spectrum 
(soHd hne) is the SDSS quasar composite from IVanden Berk et ahl (l200ll ). The quasars 
that contribute to this part of the spectrum are highly luminous high-re dshift quasa r s. Th e 
bottom spectrum (dashed line) is the HSTUV composite spectrum from lTelfer et al.l (120021 ). 
The objects contributing to this part of the spectrum are low-redshift {z < 1) quasars. The 
absence of N IV] in both composite spectra suggests that it is necessary to probe deep into 
the QLF at high-reshifts to witness the interaction of star formation with AGN activity. 



